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We present a combined measurement of the Cabibbo-Kobayashi-Maskawa matrix element \V c b\ 
and of the parameters p 2 , Ri(l), and i?2(l), which fully characterize the form factors for the 
B° — > D*~ l + vt decay in the framework of heavy-quark effective theory. The results, based on 
a selected sample of about 52,800 B° — > D*~ decays, recorded by the BABAR detector, are 
p 2 = 1.157 ± 0.094 ± 0.027, Ri{l) = 1.327 ± 0.131 ± 0.043, R 2 (l) = 0.859 ± 0.077 ± 0.021, and 
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JF(l)|F c i,| = (34.7 ± 0.4 ± 1.0) x 10~ 3 . The first error is the statistical and the second is the 
systematic uncertainty. Combining these measurements with the previous BABAR measurement of 
the form factors, which employs a different fit technique on a partial sample of the data, we improve 
the statistical precision of the result, p 2 = 1.191±0.048±0.028, Ri(l) = 1.429±0.061±0.044, R 2 (l) = 
0.827 ± 0.038 ± 0.022, and J^="(1)|V^ 6 | = (34.4 ± 0.3 ± 1.1) x 10" 3 . Using lattice calculations for the 
axial form factor .F(l), we extract \V c b\ = (37.4 ± 0.3 ± 1.2±i;l) x 10~ 3 , where the third error is 
due to the uncertainty in ^(1). We also present a measurement of the exclusive branching fraction, 
B = (4.69 ± 0.04 ± 0.34)%. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

I. INTRODUCTION 

The study of the semileptonic decay B° — > D*~ i + vi [1] 
is interesting in many respects. In the standard model, 
the rate of this weak decay is proportional to the square 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 
ment V c b, which is a measure of the weak coupling of the 
b to the c quark. This decay is also influenced by strong 
interactions. Their effect can be parameterized by two 
axial form factors A\ and A 2 , and one vector form fac- 
tor V, each of which depends on the momentum transfer 
squared q 2 of the B meson to the D* meson. The form 
of this dependence is not known a priori. In the frame- 
work of heavy-quark effective field theory (HQET) [2, 3], 
these three form factors are related to each other through 
heavy quark symmetry (HQS), but HQET leaves three 
free parameters, which must be determined by experi- 
ment. 

The extraction of \V c b\ relies on the measurement of 
differential decay rates. HQS predicts the normalization 
of the decay rate at the maximum q 2 , and \V c b\ is deter- 
mined from an extrapolation of the form factors to this 
value. The precise determination requires corrections to 
the HQS prediction for the normalization, as well as a 
measurement of the variation of the form factors near 
the maximum q 2 , where the decay rate goes to zero as 
the phase space vanishes. 

Several experiments have measured \V c b\ based on 
studies of the differential decay width for B° — > D*~ l + vi 
decays [4-10]. These analyses of the one-dimensional dif- 
ferential decay rate resulted in the measurement of only 
one of the form-factor parameters, and they relied on 
a measurement by the CLEO Collaboration [11] for the 
other two. The uncertainty in these two parameters in- 
troduces the largest systematic uncertainty in all previ- 
ous measurements of \V c b\ using this method. 

Furthermore, for measurements of \V u b\, based on both 
inclusive and exclusive B — > X u £v decays, improved 
knowledge of all form factors is important to correctly 
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describe the dominant B — > X c lv background. 

In this paper, we present a simultaneous measurement 
of \V c b\, of the branching fraction for B° — ► D*~£ + V£, 
and of the three form-factor parameters, based on mea- 
surements of three one-dimensional decay distributions. 
Thus we extend the earlier BABAR measurement [10] 
where F(1)|V^| and one of the form- factor parameters 
are measured, fully accounting for correlations between 
these one-dimensional distributions. 

We combine the results of this analysis with another 
BABAR measurement of the form factors [12], which em- 
ploys a fit to the full four-dimensional decay distribution 
on a partial sample of the data. 

The D*~ candidates are reconstructed from the D*~ 
— ► D° ir~ decays and the D° mesons are reconstructed 
in three different decay modes, K + 7t~ , i4T + 7r~7r + 7r~, and 
K + ir~ir°. Electrons or muons are paired with the D*~ 
to form signal candidates. The large data sample permits 
a precise determination of the background contributions, 
largely based on data, and thus results in smaller exper- 
imental uncertainties. 

This leads to a further reduction of the form-factor 
errors. 



II. FORMALISM 
A. Kinematic variables 




FIG. 1: Quark- level Feynman diagram for the decay B° — > 
D*~i+vt. 

The lowest-order quark-level diagram for the decay 
B° — > D*~l + vt is shown in Figure 1. This decay is 
completely characterized by four variables, namely three 
angles and the Lorentz-invariant variable w, which is lin- 
early related to q 2 and defined as 
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w 



Pb ■ Pp. 

m B m D * 



2iriBTnD* 



(1) 



where m b and m£>* are the masses of the B and the D* 
mesons (2.010 and 5.2794 GeV respectively [13]), and Pg 
and Pq* are their four-momenta. In the P rest frame 
the expression for w reduces to the Lorentz boost jd* — 

The ranges of w and q 2 are restricted by the kinematics 
of the decay, with q 2 = corresponding to 



= 1.504 



and w n 



ImBTTiD* 
1 corresponding to 

„ = {m B - m D *f = 10.69 (GeV) 2 



(2) 



(3) 



The three angular variables, shown in Figure 2, are as 
follows: 

• 9i, the angle between the direction of the lepton in 
the virtual W rest frame and the direction of the 
W in the B rest frame; 

• 6y , the angle between the direction of the D in the 
D* rest frame and the direction of the D* in the B 
rest frame; 

• x, the angle between the plane formed by the D* 
and the plane formed by the W decay. 




FIG. 2: Definition of the angles 6i, 9v, and x f° r the B° — » 
D*~£ ut decay, mediated by a vector boson W; ty s refers to 
the low momentum pion from the decay D*~ — > D°tt~ . 



B. Four-dimensional decay distribution 

The Lorentz structure of the P° — *■ D*~£ + vi decay am- 
plitude can be expressed in terms of three helicity ampli- 
tudes (P+, and Hq), which correspond to the three 
polarization states of the D*, two transverse and one 
longitudinal. For low-mass leptons, these amplitudes are 



expressed in terms of the three functions Ha^w), Rx(w), 
and R 2 (w) [2, 3]: 



Hi(w) 
where 



H T = 



m B - 



R*{l~r 2 )(w + l) 
2y/l - 2wr + r 2 



h Al {w)H % {w), (4) 



Vl-2wr + r 2 1 ± 



w 



w - 



-Pl(«/ 



Pn = 1 



1 -r 
{w - 1)(1 - R 2 {w)) 
l-r 



(5) 
(6) 



with R* = (2y'mBmfl. ) /(mB+mjj- ) and r = mp* /mj. 
The functions Ri(w) and R 2 (w) are defined in terms of 
the axial and vector form factors as, 



A 2 (w) 



R 2 (w) 



R* 2 



-Ai(w), 



E>*2 



■Ai(w). 



R*' 2 w + l 

By convention, the function Ha x (w) is defined as, 

1 2 



/lAx (w) = 



R* w + l 



Ax(w). 



(7) 
(8) 

(9) 



For w — > 1, the axial form factor A\{w) dominates, and 
in the limit of infinite 6— and c— quark masses, a single 
form factor describes the decay, the so-called Isgur-Wise 
function [14]. 

The fully differential decay rate in terms of the three 
helicity amplitudes is 

d 4 r(P° -> D*-£+v e ) _ Gm B m 2 D , 



(10) 



dwd cos 9id cos 9yd\ 8(47r) 4 
x y/w 2 - 1(1 - 2wr + r 2 )G 2 F \V cb \ 2 

x { (1- cos Of) 2 sin 2 9 v H 2 + {w) 

+ (1 + cos(9 £ ) 2 sin 2 9 v H 2 (w) 
+ 4 sin 2 9 1 cos 2 9 v H^(w) 

- 2 sin 2 9 e sin 2 9 V cos2 X H + (w)H-(w) 

— 4 sin 6^(1 — cos9g) sin#y cos#y cosx 
x H+(w)H Q (w) 

+ 4sin^(l + cos9g) sin#y cos#y cosx 

x H_{w)H (w)}. 

By integrating this decay rate over all but one of the four 
variables, w, cos9i, cos#y, or x, we obtain the four one- 
dimensional decay distributions from which we extract 
the form factors. The differential decay rate as a function 
of w is 



dT 

dw 



G 2 ? 
48^' 



<-D-+\ 



m Bo - m D ,+] Q(w)T 2 (w)\V cb \ 2 1 



(11) 
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where 

T 2 {w)g{w) = h 2 Ai (w)V^l(w + 1) 

1 



1 - 2wr + 



l + Ri(w)- 



. w 



+ 



w + 1 

and Q (w) is a known phase space factor 

1 + 4 



i + (i-Mw))^— 



g(w) = vV 2 - i(w + i) 2 



1 — 2wr 



w 



(1 -r) 2 



It is important to note that /i J 4 1 (l) = .F(l) corre- 
sponds to the Isgur-Wise function [14] at w — 1. In 
the infinite quark-mass limit, the HQS normalization 
gives = 1. Corrections to this HQS prediction 

have been calculated in the framework of lattice QCD. 
A recent calculation, performed in a quenched approx- 
imation, predicts (including a QED correction of 0.7%) 



= 0.919: 



0.035 



[15]. This value is compatible with 



estimates based on non-lattice methods [16, 17]. 



C. Form-factor parameterization 

Since HQET does not predict the functional form of 
the form factors, a parameterization is needed for their 
extraction from the data. Perfect heavy quark symmetry 
implies that R\{w) = R 2 {w) = 1, i.e., the form fac- 
tors A 2 and V are identical for all values of w and differ 
from Ai only by a simple kinematic factor. Corrections 
to this approximation have been calculated in powers of 
(^1qcd/™&) and the strong coupling constant a s . Var- 
ious parameterizations in powers of (w — 1) have been 
proposed. Among the different predictions relating the 
coefficients of the higher order terms to the linear term, 
we adopt the following expressions derived by Caprini, 
Lellouch, and Neubert [18], 



(12) 



h Al (w) = h Al (l)[l-8p 2 z + (53p 2 -I5)z 2 

-(231/? 2 - 91)z 3 ], 
R^w) = i?i(l)-0.12(w-l) + 0.05(w;-l) 2 ,(13) 
R 2 {w) = i? 2 (l) + 0.11(w7-l)-0.06(w;-l) 2 ,(14) 



where z = [Vw + 1 - y/2]/[Vw + 1 + V%\- 

The three parameters p 2 , i?i(l), and i?2(l), cannot be 
calculated; they must be extracted from data. 



III. DATA SAMPLE, RECONSTRUCTION, AND 
SIMULATION 

The data used in this analysis were recorded by the 
BaBAR detector, operating at the PEP-II asymmetric- 
energy e + e~ collider. The data sample corresponds to a 
luminosity of 79 fb _1 recorded on the T(45) resonance 
(on-resonance sample), and 9.6 flV 1 recorded at a center- 
of-mass energy 40 MeV lower (off-resonance sample). 



The BABAR detector and event reconstruction have 
been described in detail elsewhere [19, 20]. The momenta 
of charged particles are measured by a tracking system 
consisting of a five-layer silicon vertex tracker (SVT) and 
a 40-laycr drift chamber (DCH) . Charged particles of dif- 
ferent masses are distinguished by their energy loss in the 
tracking devices and by a ring-imaging Cerenkov detec- 
tor (DIRC). Electromagnetic showers from electrons and 
photons are measured in a CsI(Tl) calorimeter (EMC). 
These detector components are embedded in a 1.5-T 
magnetic field of the solenoid. Electron candidates are 
selected on the basis of the ratio of the energy detected 
in the calorimeter to the track momentum, the calorime- 
ter shower shape, the energy loss in the drift chamber, 
and the angle of the photons reconstructed in the DIRC. 
Muons are identified in a set of resistive plate chambers 
inserted in the steel flux- return of the magnet (IFR). 
Information from the IFR is combined with the track 
momentum measurement and energy deposition in the 
EMC to improve the separation of muons from charged 
hadrons. 

The electron and muon identification efficiencies and 
the probabilities to misidentify a pion, kaon, or proton as 
an electron or muon have been measured as a function of 
the laboratory momentum and the angles with clean sam- 
ples of tracks selected from data [21]. Within the accep- 
tance of the calorimeter, defined by the polar angle in the 
laboratory frame, —0.72 < cos6>i a b < 0.92, and above 1.0 
GeV, the average electron efficiency is 91%, largely inde- 
pendent of the electron momentum. The average hadron 
misidentification rate is less than 0.2%. The muon detec- 
tion extends to polar angles of —0.91 < cos#i a b < 0.95. 
For a hadron misidentification rate of typically 2.0%, the 
average efficiency is close to 65%. 

The criteria for distinguishing charged kaons from 
charged pions are chosen to maximize the efficiency while 
controlling the background, and thus differ for the decay 
modes under study. Consequently, the efficiency varies 
from 87% to 97%. The uncertainties are typically 1%. 

We determine the tracking efficiency for high- 
momentum tracks by comparing the independent in- 
formation from SVT and DCH. We compute the effi- 
ciency for low-momentum tracks reconstructed in the 
SVT alone from the angular distribution of the "slow" 
pion, 7r~, in the D*~ rest frame. We use a large sam- 
ple of D*~ — > D°tt~ , D° — » K + tt~ decays selected from 
hadronic B decays. For fixed values of the D*~ momen- 
tum, we compare the observed angular distribution to 
the one expected for the decay of a vector meson to two 
pseudoscalar mesons. This study is performed in several 
bins of the polar angle. We define the relative efficiency 
as the ratio of the observed to the expected distribution 
and parameterize its dependence on the laboratory mo- 
mentum of the ir~ . Below 100 MeV, this efficiency drops 
very steeply, and reaches zero at about 60 MeV. 

Neutral pions are reconstructed from pairs of photon 
candidates of more than 30 MeV detected in the EMC 
and assumed to originate from the interaction point. For 



7 



photon pairs with an invariant mass within 15.75 MeV 
of the nominal 7r° mass, we perform a kinematic fit, con- 
straining the mass. We require that the probability of 
the fit exceeds 1%. The efficiencies, including the EMC 
acceptance, vary between 55% and 65% for tt° energies 
ranging from 0.3 to 2.5 GeV, for a mass resolution be- 
tween 5.5 and 7.5 MeV. 

We use Monte Carlo (MC) simulation of the produc- 
tion and decay of B mesons at the T(45*) resonance and 
of the detector response [22] to estimate signal and back- 
ground efficiencies, and to extract the observed signal 
and background distributions to fit the data. We as- 
sume that the T(45) decays exclusively to BB pairs. The 
simulated sample of generic BB events corresponds to 
roughly three times the BB data sample. 

Information from studies of selected control data sam- 
ples on efficiencies and resolutions is used to improve the 
accuracy of the simulation. Comparisons of data with 
the MC simulations have revealed small differences in the 
tracking and particle detection efficiencies, which have 
been corrected for. The MC simulations include radia- 
tive effects such as bremsstrahlung in the detector mate- 
rial. QED final state radiation are modeled by PHOTOS 
[23], and decays with radiative photons are are included 
in the signal sample. 

In the MC simulations the branching fractions for 
hadronic B and D decays are based on average values 
reported in the Review of Particle Physics [13]. B° — > 
D*~£ + V£ signal events are generated with the HQET- 
bascd form factors, using the specific parameterization by 
Caprini, Lcllouch and Neubert [18]. Values of the form- 
factor parameters are taken from measurements by the 
CLEO Collaboration [11]. B — > D**iv decays, involving 
orbitally excited charm mesons, are generated according 
to the ISGW2 model [24], and decays to non- resonant 
charm states are generated following the prescription of 
Goity and Roberts [25]. 



invariant mass of the hadrons to be compatible with 
the D° mass within ±2.5 times the experimental reso- 
lution, corresponding to ±34 MeV for D° — ► K + ir~ir a 
decays and ±17 MeV for the other decays. For the de- 
cay D° — ► K + n~iT , we accept only candidates from 
portions of the Dalitz plot where the square of the de- 
cay amplitude [26] exceeds 10% of the maximum. We 
select D*~ candidates with a momentum in the range 
0.5 < pd' < 2.5 GeV. For the irj from the D*~ decay, 
the momentum in the laboratory frame must be less than 
450 MeV, and the momentum transverse to the beam 
must be greater than 50 MeV. Finally the lepton, the 
ir~ , and the D° are fitted to a common vertex using a 
constraint from the beam-beam interaction point. The 
probability for this fit is required to exceed 1%. 

In semileptonic decays, the presence of an undetected 
neutrino complicates the separation of the signal from 
background. For a signal B° decay, the D*~ and the 
charged lepton originate from the B° and the only miss- 
ing particle is a massless neutrino. The absolute value of 
the B momentum, ps, is known from the total energy in 
the event and its direction is constrained to lie on a cone 
centered on the D*~£ + momentum vector. The opening 
angle of this cone, 0b,d*i, is computed for each event, 

2E B E D ,t -m%- m 2 D , t 

cos9 b ,d»i = s — ■ (15) 

2p B Vd-i 

Here to, E, and p refer to the mass, the energy, and 
the absolute value of the momentum. The condition 
| cos6b,d*c\ < 1-0 should be fulfilled in a perfectly recon- 
structed decay. The value of w depends on the azimuthal 
angle of the B° direction, which cannot be determined. 
We therefore approximate w by the average of the four 
values of uu corresponding to the azimuthal angles 0, tt/2, 
7r, and 37r/2, as was done in [12]. This approximation re- 
sults in an average resolution for w of 0.04. 



IV. EVENT SELECTION AND BACKGROUND 
SUBTRACTION 

A. Event selection 

The reconstruction of the events and the selection of 
candidate B° — > D*~l + V£ decays are largely common to 
the earlier BABAR analysis [10] . 

We select events that contain a D*~ candidate and 
an oppositely charged electron or muon with momen- 
tum in the range 1.2 < pi < 2.4 GeV. Unless ex- 
plicitly stated otherwise, momenta are measured in the 
T(45) rest frame, which is boosted relative to the lab- 
oratory frame, /3j = 0.56. We reconstruct D*~ in 
the decay channel D*~ — > D°ir~ , with the D° decay- 
ing to K + ir~ , K + -K~n + ir~ , or K + Tr~Tr°. The tracks 
of the charged hadrons from the D° candidate are fit- 
ted to a common vertex and the candidate is rejected 
if the fit probability is less than 0.1%. We require the 



B. Background subtraction 

The background subtraction is performed separately 
for each of the four kinematic variables (w, cos 6^, cos6>y, 
and x) ! t° be used for the extraction of the form- factor 
parameters and \V c b\. For each variable, we divide the 
events into ten subsets, each corresponding to one of ten 
bins of the distribution, and the fits are performed sepa- 
rately for the ten subsets. 

The selected events are divided into six signal samples 
by separating decays into electrons and muons and the 
three D° decay modes. In addition to signal events, each 
subsample contains background events from six different 
sources: 

• combinatorial background (events from BB and 
continuum qq production in which at least one of 
the hadrons assigned to the D*~ does not originate 
from the D*~ decay); 
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continuum background (D* £ + combinations from 



• fake lepton background^ true D* combined with 
a hadron misidentified as a lepton); 

• uncorrelated background {£ + and D*~ originating 
from the decay of two different B mesons); 

• B background involving higher mass charm states, 
either B+ -> D*X£ + u £ decays (via B + -> 

D l + i>i or non-resonant B + — > D*~ir + £ + V£ 
charm states), or B° — > D*~ X£ + v^ decays, (via 



B° — > D** £ + is e or non-resonant B° 
charm states); 



D* 



n 



£+ 



• correlated background events due to the processes 
B" -> D*-t+v t with r+ -> £+X, and 5° -> 
D*-X c with V c -> £+Y. 

Except for the combinatorial background, all back- 
ground sources contain a true D*~ — ► D^tiJ decay 
and thus are expected to exhibit a peak in the Am = 
m D ,- — mjjo distribution, where m D »- and m^ are the 
reconstructed masses of the D*~ and D° candidates. 
We determine the background distributions from data, 
except for the correlated background, which amounts to 
less than 1.8% of the total sample of selected candidates. 

We determine the signal and background composition 
in two steps. First, we estimate the combinatorial, the 
continuum, and the fake lepton background from fits to 
the Am distributions (Figure 3). Second, we fix the back- 
ground levels for these three sources and determine the 
uncorrelated background and the B — > D*X£vi back- 
ground from fits to the cos6b,d*i distributions (Fig- 
ure 4). The shape and normalization of the small cor- 
related background is fixed to the MC predictions and 
based on measured branching fractions [13]. 



1. Fits to Am distributions 

To estimate the shape and normalization of the com- 
binatorial, the continuum, and fake lepton backgrounds 
from the measured Am distributions, we use in addition 
to the on-resonance data, off-resonance data, as well as a 
set of on-resonance events in which no lepton is identified 
and a charged hadron is selected to take its place. We 
refer to this data sample as the fake-lepton sample. For 
each of these data sets, we perform an unbinned maxi- 
mum likelihood fit to the Am distributions. The data are 
fitted to a sum of a peak, due to correctly reconstructed 



D* 



D 7r s decays, and a combinatorial background. 



The peak is described as a sum of three Gaussian res- 
olution functions with three mean values, three different 
widths, and two parameters that specify the contribu- 
tions of continuum and fake lepton backgrounds relative 
to the total number of events above the combinatorial 




140 



150 

Am (MeV/c 2 ) 



160 



FIG. 3: (color online) Am distributions for the decays D*~ — » 
D°nJ with D° — » K + tv~ for events in which the 7rJ track is 
reconstructed in the SVT only, a) for on-resonance, b) for off- 
resonance data in which an electron or muon was identified, 
and c) for on-resonance data in which no charged lepton was 
found in the event. The data (points with statistical errors), 
integrated over the full w range, are compared to the result 
of the simultaneous fit (solid line) to all three distributions. 
The shaded area shows the fitted combinatorial background, 
the remainder (white area) represents the sum of the signal 
and peaking backgrounds. 



background. The combinatorial background is described, 
as in [10], by an empirical function, 



F comb (Am) = — 



1 



( Am-Amo ^l / Am \ ^ 



\Arnp 



1 , (16) 



where N is the normalization, Amo refers to the kine- 
matic threshold equal to the pion mass, and ci and C2 
are free parameters. 

Since the Am resolutions and background yields de- 
pend on the D° decay mode, and on whether the low- 
momentum pion track is reconstructed in the SVT alone 
or in both SVT and DCH, the parameters describing 
the peak contributions are determined separately for the 
three subsamples corresponding to the D° decay modes, 
each divided into two classes of events, depending on the 
detection of the slow pion. 

The off-resonance distributions are scaled to the on- 
resonance luminosity, and the simulated lepton signal 
and fake lepton samples are adjusted to reproduce the de- 
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tcction efficiencies and misidcntification probabilities de- 
termined from independently selected data control sam- 
ples. For the on- and off-resonance data and for the 
fake-lepton sample, the parameters describing the shape 
and normalization of the combinatorial background dif- 
fer, and they are therefore determined separately, while 
the mean and the widths of the Gaussian functions arc 
common. 

The fits extend in Am from 138 to 165 MeV. Given 
the very large number of parameters and the many data 
subsamples and various subsets, the fits to the Am distri- 
butions are performed in several steps. Initially, the pa- 
rameters describing the combinatorial background con- 
tributions are fixed to values determined from binned x 2 
fits to simulated distributions. With these starting values 
for the combinatorial background shapes, an unbinned 
maximum likelihood fit is performed to all subsets of the 
data, both on- and off-resonance, to determine the eight 
peak-shape variables and three additional relative peak 
yields for each data subset. To improve the agreement 
with the data, the parameters describing the combinato- 
rial background distributions are then refitted, together 
with the three peak yield parameters, with the remain- 
ing parameters describing the components of the peaking 
signal and background fixed to the results of the previous 
fit. 

As an illustration for these background fits, Figure 3 
shows the Am distributions for the decays D*~ — > D°tt~ 
with D° — > K + ir~ for events in which the tt~ track is 
reconstructed in the SVT only. The data are compared 
to the results of the combined fit to three distributions, 
for selected samples of on-resonance data, off-resonance 
data, and on-resonance data without an identified lepton. 

The fitted fractions of combinatorial, fakc-lcpton, and 
continuum background events are determined for the 
peak regions in Am, which are defined as 144 to 147 
MeV for events with the slow charged pion detected in 
the SVT and DCH, and 143 to 148 MeV for decays with 
the 7t7 detected in the SVT alone. 



2. Fits to cos 0b, d* i distributions 

In a second step, a binned x 2 fit is performed to the 
cosfen«< distributions in the range -10 < COS 6b. D* I < 
5, to determine the signal contribution and the normal- 
ization of the uncorrelated and B — > D*Xlvi background 
events. Neglecting resolution and radiative effects, signal 
events meet the constraint | cos0b,d**| < 1, while the 
distribution of B — > D*Xlvt events extends below — 1, 
and the uncorrelated background events are spread over 
the entire range considered. Because of final state radi- 
ation, the events in the electron samples also contribute 
to lower values of cos0B,D*e- 

This fit is performed separately for the six signal sam- 
ples, corresponding to the D° mesons reconstructed in 
three different decay modes, divided into events with elec- 
trons or muons. Each signal sample is further divided 



into ten subsets corresponding to ten bins in one of the 
four kinematic variables. We perform the fits separately 
for each bin, with the individual shapes for the signal 
and for each of the six background sources taken from 
MC simulation. The fraction of the combinatorial, fake- 
lepton, and continuum events are taken from the Am fits 
and fixed. 

To reduce the sensitivity to statistical fluctuations, we 
require that the ratio of B — > D*Xlvg background and 
of uncorrelated background to the signal be the same 
for all three D° decay modes, with either electrons or 
muons. The cos9b,d*i distributions for the six signal 
samples and the results of the fits are shown in Figure 4. 
In total, there are 68,840 decay candidates in the range 
| cos 9b.d*i\ < 1-2. The number of selected events and 
the fractions of background events are given in Table I. 
As expected, the fake rate is about a factor ten higher for 
decays to muons than to electrons. Except for the com- 
binatorial backgrounds, which vary significantly for the 
three D° decays modes, most of the other background 
fractions are similar for the six data samples. 



V. EXTRACTION OF \V cb \ AND 
FORM-FACTOR PARAMETERS 

We determine .F(l)|Vy and the three form-factor pa- 
rameters by extending the onc-dimcnsional least-squares 
fit to the w distribution used previously [10] to a com- 
bined fit of three one-dimensional binned distributions, 
with bin-by-bin background subtraction. We have cho- 
sen this approach to avoid the statistical limitations of a 
fit to a binned four-dimensional decay distribution. 

In principle, any kinematic observable that is sensitive 
to the form-factor parameters can be used. We have ex- 
amined the sensitivity of the four kinematic variables, w, 
cos 9 1, cos9 v and \ (see Sec. II), to the form-parameters 
and found that the \ distribution is practically insen- 
sitive, and thus we select the remaining three variables. 
As done for the background estimate, we again divide the 
distributions into ten bins with equal bin size for w and 
cosfV, and varying bin size for cos 9g in order to have a 
more similar population in all the bins for this variable. 

We perform a least-squares fit to these three projected 
one-dimensional distributions to extract the form factors 
and .F(l)|T4fc|, using the independently-determined back- 
ground estimates. We account for the correlations by 
noting that the statistical covariance between the content 
of two bins in two different one-dimensional distributions 
is determined by the common number of events in these 
bins, while it is zero for bins in the same distribution. 
Since each of the three distributions that are included in 
the fit contain the same events and thus have the same 
normalization, we reduce the total number of bins used 
in the fit by two, from 30 to 28. The choice of the bins 
that are left out is arbitrary, and it has been verified that 
the fit result does not depend on this choice. 
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□ D*?v<, □ Uncorrelated D*j? 

■ D*Xj?y f ■ Correlated D*# 

■ Fake Lepton □ Continuum ■ Combinatorial 




C0S ^B,D*j? C0S ^B,D*^ 



FIG. 4: The cos9B,D*e distributions (points with statistical errors) for selected events in the six subsamples (integrated over 
the full range in w) a,b) D — > Ktt, c,d) D — > Ktt-ktt, and e,f) D — > Ktttt , compared to the result of the fits to signal and 
background contributions. 



TABLE I: Number of selected candidates and background fractions, separately for the subsamples identified by the D° decay 
mode and the charged lepton. 



Subsamples 


Ktt e 


Ktt fj, 


Ktttttt e 


K-KTT1T 


Kmv u e 


Ktytt" /x 


Number of Selected Candidates 


15144 


12083 


10259 


7831 


13224 


10299 


Signal Fraction [%] 


85.27 ±0.31 


83.27 ± 0.37 69.25 ±0.54 


66.30 ± 0.66 


75.72 ± 0.43 


74.05 ± 0.50 


Background Source 






Background 


Fractions [%] 






B -> D'Xlvi 


3.89 ±0.16 


3.93 ±0.18 


3.76 ±0.19 


3.74 ±0.22 


4.05 ±0.17 


4.08 ± 0.20 


Fake Leptons 


0.23 ±0.04 


2.45 ±0.14 


0.12 ±0.03 


2.29 ±0.17 


0.22 ± 0.04 


2.41 ±0.15 


Uncorrelated 


3.59 ±0.15 


3.11 ± 0.16 


2.71 ± 0.16 


2.27 ±0.17 


3.38 ±0.16 


3.29 ±0.18 


Correlated 


0.37 ±0.05 


0.56 ±0.07 


0.36 ± 0.06 


0.48 ±0.08 


0.35 ±0.05 


0.47 ±0.07 


Continuum 


1.80 ±0.11 


2.08 ±0.13 


1.82 ±0.13 


1.43 ±0.14 


1.83 ±0.12 


1.60 ±0.12 


Combinatorial 


4.85 ±0.18 


4.60 ± 0.20 


21.98 ±0.46 


23.49 ±0.55 


14.45 ± 0.33 


14.10 ± 0.37 



A. The least-squares fit 

The concept for this fit is an extension of the one 
introduced in the previous BABAR analysis [10]. For a 
given bin with index i, _/V l data is the total number of ob- 



served events, N i s is the estimated number of back- 
ground events, and refers to the number of MC- 
simulated signal events. The fit function can be written 
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as 



n b ,„=28n 6i „=28 , 

X2 = E I! (< ata -iVf 9 - w 2 fe 
»=i j=i ^ fe=i 



Af 



x Cr 1 (Nf ata - Nf 9 - jr W k ^j , (17) 
^ k=i ' 

where the indices i and j run over the 28 bins, and the 
index k runs over all MC-simulated events, iY 2 MC , in bin 
i. W 4 is a weight assigned to the fc-th simulated signal 
event in bin i to evaluate the expected signal yield as a 
function of the free parameters of the fit, and CV,- is the 
covariance matrix element for a pair of bins i and j . 
Each weight W k is the product of four weights, = 

wf /,fc w £ wj ,fc wf*. 

1. The factor w{*' k accounts for the dependence of 
the signal yield on the parameters to be fitted. 
The \V c b\ dependence is trivially given by the ratio 



\v cb \ 2 /\v c 



MC|2 
cb I ' 



where the denominator is the ac- 
tual value used in the simulation, derived from the 
branching fraction assumed for the decay B° — > 
D*~( + vt. Similarly, the dependence on the form- 
factor parameters, p 2 , i?i(l), and J?2(l), is given 
by the ratio of the differential decay rate (Eq. 10), 
evaluated for the values assigned during the fit and 
for the values adopted in the simulation. 

2. The factor W £ accounts for the relative normal- 
ization of data and simulated samples. It is the 
product of the following terms: 

• the ratio of the total number of BB events, 
N BB = (85.9 ± 0.9) x 10 6 and the number of 
Monte Carlo events for the final states B°B° 
and B+B-; 

• the ratios 1/[1 + /+_//oo] for B°, 
[/+_// 00 ]/[l + /+-//00] for B+, where 
/+_//oo = B{T{AS) - B+B-)/B(r(AS) - 
B°B°) = 1.037 ± 0.029 [13] is the ratio 
of number of charged to neutral B mesons 
produced; 

• the ratio of the cc to bb pair-production cross 
section; 

• the ratio of the branching fraction B(D*~ — > 
D°tt-) = 0.677 ± 0.005 [13] and of the B° life- 
time t b o = 1.530 ±0.009 ps [13] to the values 
used in the Monte Carlo simulation. 

The uncertainties on the measured quantities used 
in this weight function are not accounted for in the 
fit; their impact is studied by repeating the fit with 
their values changed by one standard deviation. 

3. The factor W^ fc is the product of the correction 
factors for efficiencies, applied on a particle by par- 
ticle basis, which accounts for the residual differ- 
ences in reconstruction and particle-identification 



efficiencies between the data and the Monte Carlo 
simulation, as a function of particle momentum and 
polar angle. 

4. The factor wf ' k accounts for potential small differ- 
ences in efficiencies among the six data subsamplcs 
and allows the adjustment of the D° branching frac- 
tions, properly accounting for their correlated sys- 
tematic uncertainties. Wf ' k is the product of sev- 
eral scale factors that are free parameters in the fit, 
each constrained to its expected value within the 
estimated experimental uncertainty. Specifically, 
to account for the uncertainty in the multiplicity- 
dependent tracking efficiency, we introduce a factor 
Wf rfe = l+N tr k{1 — 8trk), where N tr k is the number 
of charged tracks of the D*£ candidates and 5 tr k is 
constrained to 1.0 within the estimated uncertainty 
of atrk = 0.8% in the single-track efficiency. Sim- 
ilarly, multiplicative correction factors Sg, S K and 
5 n o arc introduced to adjust the efficiencies for lep- 
tons (S e for e or <5 M for fj, ), kaons, and ir° mesons, 
each within their estimated uncertainties, a t (a e for 
e ± or <jfj, for p^), <jk, and cr n . Likewise, 5b(Kit)> 
^B(A"-n-7rir) j ^(Ktt-k ) are introduced to adjust the in- 
dividual D° branching fractions within their cur- 
rent measurement uncertainties [13]. Correlations 
between the branching fraction measurements are 
taken into account by the covariance matrix Cg. 

rS,k 



As a result, W, ' can be expressed as 



The corrections to the kaon and ir° efficiency, the 
decay multiplicity, N tr k, and the D° branching 
fraction depend on the particular event k in bin 
i. 

The complete ansatz for the X 2 function used in the 
fit is 



,fc, 



i,k ci,k 



n bin =28 n bin =28 

* 2 = E E 

i=i j=i 



Nf 



k=l 



N data _ N bkg ^ y^k 

N MC 
3 

x C-/( Nf ata - Nf 9 - W j 
^ fc=i 

(l-^) 2 (l-^) 2 {l-hrk? , (f-^o) 2 

T T 1 T + 



u 



K 



cr: 



trk 



EE* 

m— 1 n=l 



B(m) X C B(mn) 5 B{n) . 



(18) 



The indices n, m refer to the three D° decay modes. The 
addition of these extra terms allows us to fit all subsam- 
ples simultaneously, while taking into account the corre- 
lated systematic uncertainties and effectively propagat- 
ing these uncertainties, via the weights, to the uncertain- 
ties on the free parameters of the fit. 
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The fit procedure has been tested on a large variety of 
Monte Carlo generated signal samples. The fits are per- 
formed for multiple samples, comparable is size to the 
data. The resulting pulls are consistent with a Gaus- 
sian distribution, with no evidence for systematic biases. 
The width of the pull distribution has been found to be 
consistent with 1, giving confidence in the uncertainty 
extracted from the fit. 



B. The covariance matrix 

A direct consequence of this method is that the covari- 
ance matrix for the measurements is not diagonal, since 
the events in different bins are not all statistically inde- 
pendent. The total covariance matrix is the sum of three 
separate matrices: one for the measured data yields, and 
one each for the estimated signal and background yields. 
The diagonal elements of the matrices are the uncertain- 
ties of the bin contents. The covariance of bins belong- 
ing to the same distribution is zero, and the covariance 
of bins from different distributions is the variance of the 
number of events that is common for the two bins. 

For the data, the covariance matrix is determined un- 
der the assumption that the data obey Poisson statistics, 
and therefore the variance of the bin or of the intersection 
of two bins is simply the number of events. The estimated 
signal matrix is built in an analogous way, where the vari- 
ance of a number of weighted events n is approximated 
by the sum of the squares of the weights, V^ i=1 n w i- 

The calculation of the background covariance matrix 
is less straightforward. The diagonal elements arc sim- 
ply the estimated variances of the measured background, 
according to the procedure described in Sec. IV. How- 
ever, the background extraction procedure does not di- 
rectly determine the number of common events in two 
bins because this procedure is based on a rather complex 
sequence of fits to shapes and event yields. The solu- 
tion adopted here is to use the number of common back- 
ground events in two bins as predicted by the simulation, 
corrected for tracking and particle identification (PID) ef- 
ficiencies and adjustments to account for the background 
estimates. This is done in such a way that for each bin of 
one of the kinematic observables the background is equal 
to the data-based background estimate. The choice of 
the observable (w, cos cos 9y , or x) is arbitrary, and 
the variation of the results with this choice is used to 
evaluate the systematic uncertainty introduced by this 
procedure. 

Since the total number of background events estimated 
for the four different distributions is not exactly the same, 
we average the results obtained for the four distributions. 
The spread of the background normalization values is 
found to be almost twice as large as the estimated uncer- 
tainty from the error propagation. This can be explained 
by the fact that the uncertainties of the assumed signal 
and background shapes in the Am fits are not accounted 
for in this propagation. We derive an estimate of this 



additional uncertainty using the number of background 
events determined from the fits to the background dis- 
tributions. They are very similar for w and cos^, and 
for cos Qy and \- The estimate is given by the average 
of the minimum and maximum difference between these 
two groups. 



VI. THE FIT RESULTS AND SYSTEMATIC 
UNCERTAINTIES 

A. Results of the fit 

The results of the simultaneous fit to the three one- 
dimensional distributions are presented in Table II. The 
stated uncertainties on p 2 , Ri(l), R 2 (l) and JF(1)|V C 6| 
are taken from the MINUIT minimization program [27]. 
Among the form-factor parameters, the correlations are 
quite large, but their correlation with ^(1)114^1 is less 
than 0.23. All the S parameters are found to be consis- 
tent with their nominal value of one within their uncer- 
tainties. 

Figure 5 compares the one-dimensional projections for 
the kinematic observables with the results of the fit, de- 
tailing the signal and background contributions. To pro- 
vide an additional check on the background estimation, 
we also show the x projection. The goodness of the fit, 
ignoring that the measurements used are not all indepen- 
dent, can be stated as x 2 /d.o.f. = 23.8/24, corresponding 
to a probability of 47.3%. 

TABLE II: The results of the fit: The parameters, their un- 
certainties, and the off-diagonal elements of the correlation 
matrix. The stated uncertainties include systematic ones that 
are not common to all events and are therefore included in the 
fit (see Sect. VIB1). 



Fitted Correlations 
Parameters Values p 2 ^i(l) ^2(1) 
-^(1)1 V cb | (34.67 ± 0.86) x 10" 3 -0.001 -0.196 +0.141 
p 2 1.157 ±0.095 +0.867 -0.924 
Ri(l) 1.327 + 0.131 -0.928 
R 2 (l) 0.859 ± 0.077 



The B° — > D* l + vi branching fraction, obtained by 
integrating the data over the full phase space, is found 
to be B(B° -» D*~£+u e ) = (4.72 ± 0.05)%. 

As a cross check, we perform the fit separately for the 
six subsamples. The quality of the fits is generally good, 
and the results, shown in Table III, agree within the un- 
certainties obtained, with the possible exception of the 
Kmnr e sample. Detailed investigations of the back- 
ground estimates and fits for this sample did not reveal 
any anomalies. 
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TABLE III: Results of fits performed separately for the six subsamples corresponding to each combination of three D° decay 
modes and the charged lepton. The uncertainties represent the total uncertainty of the fit, except for where it is 

split into the statistical and the systematic contribution included in the fit. 



Subsample 



R 2 (l) 



■F(l)|Vct| x 10* x /d-o.f. 



Kire 0.971 ±0.163 1.166 ±0.182 0.977 ±0.107 34.76 ±0.61 ±0.61 23.9/24 

Kivn 1.013 ±0.175 1.193 ±0.206 0.922 ±0.123 34.55 ±0.66 ±0.65 37.9/24 

Ktyttty e 1.581 ± 0.151 2.043 ± 0.384 0.405 ± 0.232 33.30 ± 1.27 ± 0.96 15.6/24 

KTv-KTVfi 1. 146 ±0.258 1.156 ±0.351 0.946 ±0.197 34.14 ± 1.10 ±0.98 28.0/24 

Kit-it e 1.042 ± 0.165 1.217 ± 0.206 0.926 ± 0.118 34.86 ± 0.64 ± 1.46 26.9/24 

/i 1.170 ±0.155 1.439 ±0.228 0.838 ±0.131 34.38 ± 0.74 ± 1.46 24.8/24 
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FIG. 5: Comparison of the measured distributions (data points) a) w, b) cosdt, c) cos#y, and d) x, with the result of the fit, 
shown as the sum of the fitted signal yield and the estimated background distributions. The statistical uncertainties of the 
data are too small to be visible. 



B. Systematic uncertainties 



Uncertainties included in the fit 



A summary of statistical and systematic uncertainties 
on the measured parameters is presented in Table IV, 
including the breakdown of those for the measurement of 
the B° — > D*~ l + vi branching fraction. 



The uncertainty of the parameters resulting from the 
fit is not purely statistical, since the systematic uncer- 
tainty sources that are not common to all events are ac- 
counted for in the fit through the 8 parameters in the 
weights Wf' fc . As described above, these weights account 
for residual uncertainties in the lepton and hadron iden- 
tification, the charged particle tracking and 7r° efficien- 
cies, and the individual D° branching fractions. We can 
determine the statistical uncertainties of the fit by re- 
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TABLE IV: Breakdown of statistical and systematic uncertainties. 





„2 
P 


rtl(l) 


R»{1) T{l)\V cb \ x 10 3 - 


-» U i: T //|J X 11) 


Statistical Error 


0.094 


0.131 


0.077 


0.41 


0.05 


nTr\ j_ i ■ k>/ r~»U\ 

FID, tracking, B(D ) 


0.003 


0.006 


0.002 


0.75 


0.21 


Soft-7r efficiency 


0.013 


0.005 


0.001 


0.46 


0.18 


D*l vertex fit 


0.014 


0.010 


0.008 


0.06 


0.06 


I?-momentum variation 


U.Ulo 


U.U4U 


0.017 


0.29 


n 1 a 
U.14 


Radiative corrections 


u.uuo 


fl nfi/1 
U.UU4 


0.000 


0.19 


U.U / 


D** composition 


0.011 


0.008 


0.009 


0.10 


0.07 


Background estimates 


0.006 


0.004 


0.002 


0.04 


0.04 


Partial Systematic Error 


0.027 


0.043 


0.021 


0.95 


0.33 


B u lifetime 








0.10 


0.03 


BB normalization 








0.19 


0.05 


B(D* D°it) 








0.13 


0.04 


/+-//00 


0.003 


0.003 


0.002 


0.25 


0.07 


Total Systematic Error 


0.027 0.043 0.021 


1.01 


0.34 



pcating the fit with all <5 parameters fixed at their fitted 
value. The systematic uncertainties are then obtained 
by subtracting the statistical covariance matrix from the 
total covariance matrix. While the uncertainties related 
to the detector performance are relatively small for the 
form-factor parameters, they are dominant for JF(l)|V c b| 
(2.5%) and the branching fraction (5.2%). 



2. Soft pion efficiency 

A major source of uncertainty on .F(l)|T4f,| is the re- 
construction efficiency for the low-momentum pion from 
the D*~ decay, since it is highly correlated with the D*~ 
momentum and thereby with w. The functions param- 
eterizing the efficiency for data and MC simulation arc 
consistent within the statistical uncertainties. To assess 
the systematic uncertainty on | V c b | , we vary the parame- 
ters of the efficiency function by their uncertainty, includ- 
ing correlations. We add in quadrature the uncertainty in 
the absolute scale as determined from higher-momentum 
tracks reconstructed in both the SVT and the DCH. The 
resulting systematic uncertainty on \V c b\ is 1.3 %. 



3. D*l vertex reconstruction efficiency 

The uncertainties from the D*l vertex reconstruction 
have been evaluated by observing the impact of changes 
in the standard vertex fit procedure. First, we remove 
the constraint on the average position of the beam-beam 
interaction point, and second, we remove the lepton track 
from the vertex fit. We take the larger of the observed 
variations of the parameters as the systematic uncer- 
tainty. 



4- B momentum 

The B momentum is impacted by small changes in the 
energies of the two colliding beams. In the Monte Carlo 
simulation the beam energies are assumed to be constant. 
We have examined the impact of these variations on the 
measured distributions and have concluded that we can 
account for this difference by rescaling the cos9b,d*i val- 
ues in the simulation by a factor of 0.97. Half of the 
observed relative change of the fitted parameters coming 
from the adjustment is assumed as the systematic uncer- 
tainty due to this effect. This is the largest systematic 
uncertainty on i?i(l) and i?2(l)- 



5. Radiative corrections 

Radiative corrections to the B° — ► D*~i + vi decays 
are simulated by PHOTOS 2.0 [23], which describes the 
final state photon radiation (FSR) up to O(ot 2 ). In the 
event reconstruction no attempt is made to recover pho- 
tons emitted in the decay. The simulated prediction of 
the reconstructed kinematic variables is sensitive to the 
details of the radiative corrections. This is particularly 
important for electrons, for which FSR results in the long 
tail below —1 in the cos9b,D'1 distribution. 

At present, no detailed calculation of the full 0(a) 
radiative corrections to the B° — > D*~£ + vg decay is 
available. Recently, a new 0(a) calculation of radiative 
corrections in K° — > ir~e + is e decays has become avail- 
able [28], which allows detailed comparisons of the ra- 
diated photons with data and with PHOTOS calculations. 
These new calculations agree well with kaon data. From 
the comparison with PHOTOS, it is evident that the radi- 
ated photon energy spectrum is quite well reproduced by 
PHOTOS, but the photon emission angle with respect to 
the electron differs significantly. 

To assess the systematic uncertainty due to the imper- 
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feet treatment of the radiative corrections in B decays, 
we have used the comparison presented in Ref. [28] and 
reweighted the simulated decays to reproduce the photon 
angular distribution for photons above 10 MeV in the B 
rest frame. The effect of the reweighting has been used as 
an estimate of the systematic uncertainty. The possible 
impact of the radiative effects on the hadronic interac- 
tions, i.e. the form factors, is unknown and therefore 
not considered. 



6. B — > D* Xivi background description 

The semileptonic branching fraction and form factors 
for the four higher-mass charm states, D** , and for 
non-resonant production of Dtt and D*ir are not well 
known. The shapes of these different components of 
the B — > D*X£vt background are taken from simula- 
tion. Their relative yield is obtained from the fit to the 
cos6b,d*i distribution. To account for the uncertainty of 
the composition of this background, the fits have been re- 
peated using only one of the contributions at a time. The 
study was done only for the D° — > K + tt~ subsamplc; it 
is assumed to be valid also for the other subsamplcs. The 
estimated uncertainty on the fit parameters is taken as 
half of the biggest change observed. 

7. Background estimates 

As explained in Sec. VB the background covariance 
matrix is built using the measured background shape in 
one observable, and the simulation information for the 
others. The choice of the observable is arbitrary. The 
systematic uncertainty is evaluated by comparing results 
for the four kinematic observables. The maximum ob- 
served variation with respect to the standard fit result is 
taken as the estimate for this uncertainty. 

8. Global normalization factors 

There are several quantities that affect only the overall 
normalization of the data, and thus not the form-factor 
parameters. Their contribution to the systematic uncer- 
tainty on ,F(1)|V^&| and the branching fraction are listed 
in the bottom half of Table IV. They are: the B° lifetime 
(t b o = 1.530 ± 0.009 J13]), the D* 1 -> D°n- branching 
fraction (B(D*~ D tt~) = 67.7 ±0.5%), and N BB , 
the number of BB events in the total data sample. The 
systematic uncertainty on N BB is 1.1%. 

The effect of the uncertainties in the D° branching 
fractions has already been discussed; it is subsample spe- 
cific, but it affects all parameters because it changes the 
fraction of signal events from different D° decays. 

The uncertainty on the ratio f + _/f 00 = B(T(4S) — > 
B+B-)/B{T{AS) -» B°B°) = 1.037 ±0.029 [13] affects 
both the absolute number of measured B° — > D*~t + h>i 



decays and the ratio of background events from B° and 
B^ decays in the simulation. This second aspect influ- 
ences the cos8b,d*£ distributions and therefore the back- 
ground determination. For this reason, this uncertainty 
affects also the form-factor parameters. The systematic 
uncertainty is equated with the observed change in the 
fit parameters for a one-standard deviation change in the 
value of /+_//oo- 

VII. CONCLUSIONS 
A. Summary of results of this analysis 

A sample of about 52,800 fully reconstructed B° — ► 
D*~£ + ve decays recorded by the BABAR detector has 
been analyzed to extract both and the form- 

factor parameters, p 2 , -Ri(l) and i? 2 (l)j m the Caprini- 
Lellouch-Neubert parameterization [18]. The D*~ candi- 
dates are reconstructed from the D*~ — > D° n~ decays 
and the D° mesons are reconstructed in three different 
decay modes, K + n~ , K + tt^tt + 'k^ , and K + -k~tt°. Elec- 
trons or muons are paired with the D*~ to form signal 
candidates. The large data sample has permitted a more 
precise determination of the background contributions, 
largely based on data, and thus has resulted in smaller 
experimental uncertainties. 

The results of the simultaneous fit to three one- 
dimensional projections of the decay rate are 

•F(l)|V c6 | = (34.7 ± 0.4 ± 1.0) x 10" 3 

p 2 = 1.157 ±0.094 ±0.027 

i?i(l) = 1.327 ±0.131 ±0.043 

R 2 (l) = 0.859 ±0.077 ±0.021. 

The stated uncertainties of the measurement here are the 
statistical and the total systematic one. The simultane- 
ous fit to the three distribution reduces the uncertainty 
due to the form- factor parameters. 

Using an unquenched lattice calculation giving = 
0.919^3° [15] results in the following value for |U c6 |, 

|V c6 | = (37.8 ± 0.4 ± 1.1±1;2) x 10- 3 . 

Here the third error is due to the theoretical uncertainty 
in .F(l). Figure 6 shows the measured decay rate, inte- 
grated over angles, V^|, as well as the fitted theo- 
retical w dependence (see Eq. 11). 

The branching fraction for the decay B° — > D*~ ( + vn 

is 

B(B° -» D*-t+i/ t ) = (4.72 ± 0.05 ± 0.34)%. 



B. Combination of results with the previous BABAR 
measurement of the form-factor parameters 

The BABAR collaboration recently published a mea- 
surement [12] of the same form- factor parameters for 
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w 

FIG. 6: The measured w dependence of J r («;)|V" C 6| (data 
points) compared to the theoretical function with the fitted 
parameters (solid line). The experimental uncertainties are 
too small to be visible. 



B° — ► D*~ t + Vf decays based on an unbinned maximum- 
likelihood fit to the four-dimensional decay distribution 
(Eq. 10). This fit is sensitive to the interference of the 
three hclicity amplitudes and thus results in significant 
smaller uncertainties on the form-factor parameters. The 
fit does not attempt an absolute normalization of the de- 
cays, and thus is not sensitive to .F(l)|T4f,|. It resulted in 
p 2 = 1.145±0.066±0.035,i?i(l) = 1.396 ±0.070 ±0.027, 
and R 2 (l) = 0.885 ± 0.046 ± 0.013. 

We combine the two BABAR measurements of the form- 
factor parameters, taking into account the correlation be- 
tween them, and obtain 

■^(1)1^1 = (34.4 ± 0.3 ± 1.1) x 10" 3 

p 2 = 1.191 ±0.048 ±0.028 

Ri(l) = 1.429 ±0.061 ±0.044 

R 2 (l) = 0.827 ±0.038 ± 0.022. 

Compared to the analysis presented in this paper, the 
combined result has significantly smaller statistical un- 
certainties of the form-factor parameters. The event sam- 
ple and the sample of Monte Carlo simulated events used 
in Ref. [12] are a subset of the one used in the present 
analysis, namely about 15,000 selected B° candidates 
with D° — > K + 7r~ decays combined with electrons. Ex- 
cept for the selection of the D° decay, the event selection 
and the determination of the backgrounds shapes and the 
signal extraction are almost identical for the two analy- 
ses. Therefore, all the detector-related systematic uncer- 
tainties should be the same, as well as the uncertainties 
from the background models and input parameters like 
the branching fractions. Thus, we retain the systematic 
measurement uncertainties established in this paper. The 
combined statistical errors are still larger than the total 
systematic uncertainties, but not by a large factor. An 
upper limit for the correlation between the two measure- 
ments has been estimated on the basis of the ratio of the 
uncertainties, and is found to be less than 0.45. 

The correlation coefficients for the combined measure- 



ments are 



P(p 2 ,i?i(l)) 


= +71% 


P(p 2 ,i? 2 (l)) 


= -83% 


p{p 2 ,T{l)\V cb \) 


= +27% 


p(R 1 (l),R 2 (l)) 


= -84% 


p(i2i(l), ^(1)1^1) 


= -39% 


p(R 2 {l),T{l)\V cb \) 


= +22%. 



Figure 7 shows the correlations between the fitted vari- 
ables and their uncertainties, both for the present analy- 
sis and for the combined result with Ref. [12]. The con- 
tours correspond to A% 2 = 1, i.e. 39% CL. The cor- 
relations between the form-factor parameters are quite 
large, but their correlation with JF(l)|V" c fc| is less than 
0.4, and the sign of the coefficients differ, resulting in a 
much reduced overall dependence of .F(l)|V c b| on these 
form factors. 

Using the same lattice calculation for [15], we 

obtain an improved value for \V c b\, 

\V cb \ = (37.4 ± 0.3 ± 1.2 x 10~ 3 , 

where the third error reflects the current uncertainty on 

The corresponding branching fraction of the decay 
B° -> D*-£+v t is 

B(B° -» D*-t + v t ) = (4.69 ± 0.04 ± 0.34)%. 

The combined results of the two BABAR analyses super- 
sede all previous BABAR measurements of the form-factor 
parameters, of the exclusive branching fraction for the 
B° — ► D*~£ + vg decay, and of \V cb \ extracted from this 
decay. 

The value of the branching fraction presented here is 
smaller than the average of previous measurements [13]. 
This measurement combined with B(B° — > D^lvg) = 
(2.08 ± 0.18)% [13] represents only (65 ± 7)% of the to- 
tal branching fraction for the B° — > XJLvt decays. The 
remaining fraction of 35% is expected to involve higher- 
mass charm states. The branching fractions for decays to 
these individual higher-mass states are not well known, 
in particular those involving broad resonances or non- 
resonant LK*'**)^ states [29, 30]. 

The combination of the two BABAR measurements re- 
sults in a further reduction of the form-factor uncertain- 
ties compared to the previous BABAR analysis [12], for 
which the uncertainties on -Ri(l) and i?2(l) had already 
been reduced by a factor of four or more, compared to 
the CLEO measurement [11]. The uncertainty on p 2 has 
also been reduced, by a factor of five with respect to 
the BABAR measurement in Ref. [10]. The correlation 
between .F(l)|Vd,| and p 2 , which was sizable for all pre- 
vious measurements, has been reduced significantly, and 
this also leads to a smaller uncertainty for \V cb \. 

The resulting value of \ V cb \ is fully compatible with the 
earlier BABAR measurement [10], and most earlier mea- 
surements [13], but it is significantly smaller than the 
CLEO measurement [31]. 
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FIG. 7: (color online) Correlations between fitted variables and their uncertainties, both for the present analysis with the 
statistical uncertainties (dotted line, central value marked with an open circle) and for the combined result with the total 
experimental uncertainties, i.e. the statistical and systematic uncertainties combined (solid line, central value marked with a 
solid dot), a) p 2 -Ri{l), b) p 2 -R 2 {l), c) p 2 -T(l)\V cb \ x 10 3 , d) e) R x {l)-F{l)\V cb \ x 10 3 , and f) R 2 {l)-T(l)\V cb \ x 10 3 

projections. The contours correspond to Ax 2 = 1, i.e. 39% CL. 



\V c b\ can also be extracted from inclusive B — > X c tvi 
decays. Recent measurements of \V c b\ rely on moments 
of the electron energy and the hadron mass spectrum, 
combined with moments of the energy photon spectrum 
in inclusive B — > X s j decays. Here X c and X s refer 
to charm and strange hadronic states, resonant or non- 
resonant. Global fits to such moments, as a function of 
the minimum lepton and photon energy have been per- 
formed in terms of two different QCD calculations, one in 
the so-called IS scheme [32] and the other in the kinetic 
scheme [33] . The results are in very good agreement. The 
weighted average is \V cb \ = (41.7 ± 0.7) x 10 -3 [13]; the 
experimental and theoretical uncertainties are compara- 
ble in size. The experimental techniques and the theo- 
retical calculations for exclusive and inclusive decays are 
completely different and uncorrelated. Given the sizable 
experimental error and the large uncertainty of the form- 
factor normalization of the exclusive measurement, the 
results are consistent. This should give us confidence in 
the experimental methods employed and the theoretical 
calculations, but also an incentive to improve on both. 
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